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MSX VERSUS IRAS TWO-COLOR DIAGRAMS AND THE CSE-SEQUENCE OF OXYGEN-RICH LATE-TYPE STARS 
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ABSTRACT 

We present MSX two-color diagrams that can be used to characterize circumstellar environments of sources 
with good quality MSX colors in terms of IRAS color regions for oxygen-rich stars. With these diagrams we 
aim to provide a new tool that can be used to study circumstellar environments and to improve detection rates 
for targeted surveys for circumstellar maser emission similar to the IRAS two-color diagram. This new tool 
is especially useful for regions in the sky where IRAS was confused, in particular in the Galactic plane and 
bulge region. Unfortunately, using MSX colors alone does not allow to distinguish between carbon-rich and 
oxygen-rich objects. An application of this tool on 86 GHz SiO masers shows that for this type of masers an 
instantaneous detection rate of 60% to 80% can be achieved if target sources are selected according to MSX 
color (region). 

Our investigations may have revealed an error in the MSX point source catalog version 2.3. That is, the 
photometry of the 21.3 /im (MSX E filter) band for most weak 8.28 /im (or MSX A filter) band sources seems 
off by about a factor two (0.5-1 magnitude too bright). 

Subject headings: masers — catalogs — stars: AGB and post-AGB — circumstellar matter — stars: late-type 
— infrared: stars 



1. INTRODUCTION 

The two-color diagram (2CD) of stellar dust/gas en- 
velope sources in the Infrared Astronomical Satellite 
(IRAS) point source catalog (PSC) has been very use- 
ful in characterizing the circumstellar environment of late- 
type stars, typically thought to be at the top of the 
Asymptotic Giant Branch (AGB) (iNeugebauer etal.l l 19841 
Beich man et al.lll985tloInon et al.lll984tlHacking et aljl985t 



the sense that due to the low angular resolution, the Galac- 
tic plane (and in particular the bulge) survey region suffers 
from confusion. This is unfortunate as the majority of evolved 
stars (but also e.g. young stellar objects), are found in the 
plane. Therefore subsequent infrared surveys using the In- 
frared Space Observatory (ISO, with the ISOGAL survey), 
the Midcourse Space Experiment (MSX, with the Galactic 
Plane survey), and the Spitzer Space Telescope (with the 
GLIMPSE and MIPSGAL surveys) have concentrated on the 
Galactic plane and bulge with much better angular resolution 3 
dOmont et all 120031: iPrice eTafl l200lt iBeniamin et al.l l200l 
Carey et al. 2009^7 However, in correcting for IRAS' short- 
coming, areas of the sky which overlap with the non-confused 
IRAS observations are minimal. Also, these other surveys 
have typically used other infrared wavebands and filters than 
IRAS. When the sources of interest are the intrinsically vari- 
able late-type stars, a direct comparison of colors between 
IRAS and e.g. MSX is, in principle, very difficult if not im- 
possible. 

Nevertheless, with the assumption that colors do not differ 
too much between observing epochs, i.e. stellar phase, and 
thus that results of a comparison between subsamples in a sta- 
tistical sense can be interpreted, an attempt at a comparison 
between IRAS colors and other infrared colors can be made. 
Because the PSC of MIPSGAL and AKARI were not avail- 
able in the summer of 2008 (when we began this study), and 
because ISOGAL and GLIMPSE were limited in spatial cov- 
erage and wavelength range respectively, we here present a 
comparison of IRAS and MSX colors. This comparison pro- 
vides a tool to study, in the first instance, oxygen-rich AGB 
and RGB stars in MSX colors, similar to the IRAS 2CD. We 
note that the IRAS 2CD cannot fundamentally be replicated 
in MSX colors, because the MSX data do not extend to a long 
enough wavelength so that the vertical axis of the IRAS 2CD 
can be approximated with MSX data. Most of the analysis and 
Electronic address: |lsjouwerman@nrao.edu, stephanie.m.capen@enc.edu, mclau ssea@6tia&a>11 that follows concentrates on the horizontal extent 

'National Radio Astronomy Observatory, P.O. Box O, Socorro, NM of the IRAS 2CD. 

87801 

2 Eastern Nazarene College, 23 East Elm Avenue, Quincy, MA 02170 3 The AKARI mission h an all . sky survey jMurakami et alllDOl) . 



van der Veen & Habing 1988). The IRAS 60^m/25/im versus 
25/im/12/im two-color diagram statistically distinguishes be- 
tween carbon-rich and oxygen-rich objects. In the latter color 
it also indicates the thickness of the circumstellar shell based 
on the appearance of the 9.7 /im silicate feature in the individ- 
ual 7.7 to 22.6 um low-resolution sp ectra (LRS) (Oln on et ail 
[19861 Ivan der Veen & Habingll 19881) . 

For many studies since the IRAS mission, the IRAS 2CD 
has thus been the starting point for color-selection of samples; 
not only to select AGB or Red Giant Branch (RGB) stars, but 
also to select e.g. proto-planetary nebulae, young stellar ob- 
jects and starburst galaxies. The IRAS 2CD has been revolu- 
tionary in the study of the late evolutionary stages of evolved 
stars and chemical enrichment of the Galaxy due to the en- 
hanced mass-loss of these stars. Moreover, the 2CD selection 
is very effective in finding circumstellar masers, and hence 
line-of-sight velocities of stars that then can be used for kine- 
matic studies of the Gal axy. For example, objects in IRAS 
color region Illb and IV (Ivan der Veen & Hab ing 1988|) have 
been targeted for 1612 MHz OH maser emission, whereas ob- 
jects in region II usually have a circumstellar envelope (CSE) 
too thin to form and sustain these OH masers. On the other 
hand, objects in region II may be CSE s in which Si O maser 
emission is more likely to exist (e.g.. lHabing||1996l and the 
extensive list of references therein). 

The IRAS all-sky survey has a significant shortcoming in 
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FIG. 1. — Positional offsets of the 2543 cross-identification s between the 
IRAS and MSX catalogs within 10" radius using the van der Veen & Habing 
1 1988) IRAS 2CD boundaries and with at least one good quality MSX color. 



As an example of the capability of this comparison, we 
apply our MSX color tool to investigate the occurrence of 
86 GHz SiO masers in color-selected sources in the Galac- 
tic bulge. It is our intention to derive similar comparisons for 
Spitzer (combined GLIMPSE and MIPSGAL) and AKARI 
colors once their point source catalogs are public, and to in- 
vestigate the occurrence of other CSE masers as a function of 
color in a similar way using these new tools. 

Section [2] describes the selection of the data sets; the 
sources that are in both the IRAS and MSX PSCs but ex- 
cluding sources that could confuse or contaminate a statis- 
tical analysis. As most masers originate from oxygen-bearing 
molecules, we would like to be able to separate the carbon- 
rich from the oxygen-rich CSEs based purely on MSX col- 
ors. But before investigating such a separation in Sect. [3] 
we compare the "evolutionary curve" of oxygen-rich stars 
with circumstellar mate r ial in t he IRAS 2CD as described by 
Ivan der Veen & Habingl (119881) with MSX colors. Note that 
this "evolutionary curve" is more a circumstellar envelope 
opacity measure than a strict stellar mass-loss evolution curve. 
Therefore we will refer to this curve as the "CSE-sequence" 
for oxygen-rich late-type stars in the remainder of this paper. 
The resulting MSX 2CDs are discussed in Sect. [4] and an ex- 
ample of their scientific usefulness is shown in Sect. [5] Future 
plans are collected in Sect. |6]before closing with a summary. 

2. IRAS AND MSX DATA SELECTION 
2.1. Cross-identification between IRAS and MSX sources 

For the cross-identification between IRAS and MSX 
sources we used the positions from version 2.0 of the 
IRAS PSC (245889 sources) and Faint Source Catalog (FSC; 
173044 sources) and version 2.3 from the MSX(6C) P SC 
(440487 sources) dBeichman et al.lll988l;lEgan et al.ll2003l) . 

Because we were very wary of contamination with incorrect 
cross-identifications, but on the other hand needed a data set 
of cross-identification as large as possible, we compromised 
between mission position accuracy and search radius. That is, 
IRAS PSC 12 /im positions typically are accurate to ~ 7-9", 
within a beam (angular resolution) of ~ 1 x5', and MSX PSC 
8.28 /im positi ons typically are accura t e to ~ 2", within a 
beam of ~ 1 8" dNeugebauer et alJI 1 9841; Backing et alJI 1 985t 
lOiha et al.l20"07tlBeichman et al.ll988tlEgan et al.ll2003l) . We 
therefore used a search radius of 10"; the arithmetic mean of 
the beam and accuracy of the MSX mission, which is the more 
restrictive between IRAS and MSX. We used the same 10" 
to cross-identify the MSX PSC with the IRAS FSC for two 
reasons. First, this expands the data set with more sources 
to compare IRAS and MSX colors. Second, it also reveals 
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FIG. 2. — The "traditional" IR AS 2CD with uncorrect ed IRAS colors, color 
regions and color ranges as in lvan der Veen & Habing 1 1988, their Fig. 5b). 
The (solid) black body line, with temperature ticks, and the (dashed) "evolu- 
tionary track" are also shown. The "evolutionary track" is referred to as the 
"CSE-sequence" in the remainder of this paper. The data (points) are all 2543 
good quality IRAS sources with a likely unique MSX counterpart within 10" 
radius for which at least one good quality MSX color can be derived. 



those IRAS-MSX PSC cross-identifications that could be con- 
fused with faint IRAS sources. Our data set thus excludes the 
MSX sources that had more than one possible IRAS counter- 
part within 10", whether or not from the IRAS PSC or FSC, 
but also includes the (single) fainter IRAS counterparts. In to- 
tal we found 45895 MSX sources identified with single IRAS 
PSC sources and 620 identified with single IRAS FSC sources 
(46515 total). Multiple identifications within 10" were found 
for 482 MSX sources, whereas 393490 MSX sources had no 
counterpart in the IRAS PSC nor in the IRAS FSC within 10" 
radius. 

Two MSX sources were removed because they could be (at 
6.64" or 5.54") identified with the same FSC source. Sim- 
ilarly, we removed another 76 MSX sources for which the 
identification with a unique IRAS PSC source within 10" 
was ambiguous. To be consistent with the multiple identi- 
fications above, this means removing the MSX source even 
in the case where one possible identification was within 3" 
(2.87") and the other was as far as (almost) 10" (9.60") away. 
From this sub-total of 46437 sources which match in posi- 
tion, we removed 26743 cross-identifications for which the 
MSX sources had fewer than two reliable flux densities (i.e., 
"quality numbers" for the photometry which were not 3 or 
4), i.e., for which not a single reliable, good quality MSX 
color could be calculated. Another 16194 IRAS sources 
were removed for which any of the 12, 25 or 60 /itm fluxes 
were unreliable (quality not 3) as in that case we could not 
place them in the IRAS 2CD diagram reliably. Finally, an- 
other 957 sources were removed as th ey fall outside the tradi- 
tional |vindir^en & Habing (1981) IRAS 2CD color mag- 
nitude range limits (their Fig. 5a, 5b: -2 < [12] -[25] < 2 
and -3 < [25] -[60] < 1; see next subsection). The re- 
maini ng catalog of good quality IR AS sources in the tradi- 
tional lvan der Veen & Hab ing ( 1988) IRAS 2CD with a likely 
unique MSX counterpart within 10" radius for which at least 
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Fig. 3. — Top: uncorrected IRAS [12]— [25] color versus combi- 
nations of the uncorrected good quality MSX colors for the 2455 
sources (2543 less 88 from IRAS region VIII) in Figs.[TJ[2] For the 
MSX [A]-[D] color, the [A]-[E] color and the [C]-[E] color there 
is a good correspondence to the IRAS [12]— [25] color, with [A]-[E] 
having a larger color range (a bigger stretch, steepest slope: TableQJ. 
Bottom: same for the IRAS [25]-[60] color, where no one-to-one 
correspondence between the IRAS [25]-[60] color can be deduced 
for any MSX color. 



one good quality MSX color could be derived is 2543. Only 
two sources are from the FSC. Figure Q] shows the histogram 
of the difference in position between the final 2543 IRAS and 
MSX associations. From the peaked distribution around 2.5" 
(with ~ 18" angular resolution) we conclude that the associa- 
tions are very likely proper cross-identifications. The location 
of these remaining cross-identified IRAS/MSX sources in the 
IRAS 2CD diagram is shown in Fig. [2] 

2.2. Flux density and color definitions 

To be consistent with the original IRAS 2CD nomencla- 
ture we follow the (uncorrected ) color magnitude definition 
used by lvan der Veen & Habingl (11988). i.e. color magnitude 
[a] - [f3] = —2.5 x log 10 (F Q ,/F ( g) where F x is the measured 
flux density in Jansky (Jy) at waveband \ taken from the 
IR AS PSC/FSC version 2.0 or from the M SX PSC version 
2.3 dBeichman et alj!988tlEgan et al]2003l) . For IRAS wave- 
bands we will use the numbers 12, 25 and 60 to designate 
the waveband; the center wavelength in /im of the filter pass- 
band. For MSX we will use the capital letters A, C, D and 
E for the wavebands centered at 8.28, 12.13, 14.65 and 21.3 
/im respectively. If we call a source "blue" or "bluer" we 
refer to the smaller or more negative value(s) of a color; 
relatively larger or more positive color values are "red" or 
"re dder". We also use the IRAS color regions determined 
by Ivan der Veen & Habingl (fl988). That is, statistically, the 
oxygen-rich stars are located in the IRAS color regions iden- 
tified with capital Roman numerals I, II, Ilia, Illb, IV and 
V along the CSE-sequence, and carbon-rich stars typically 
are located in the IRAS color regions Via, VIb and VII 





i | i i i i | i 




1 : 




















i 1 1 1 1 i 1 1 1 1 1 1 1 i 












1 1 1 1 M 1 ) 1 1 1 1 1 1 1 






i 'I'^jiftfH^'^^ v ' : ' ; '^ 










,],,,,[, 







IRAS [12]-[25] 



LOG (IRAS 12) 



FIG. 4.— IRAS [12]-[25] color (left) and IRAS 12 fj,m flux (right) ver- 
sus the difference between IRAS [12]-[25] color and MSX [A]-[D] color 
(top), MSX [A]-[E] color (middle) and MSX [C]-[E] color (bottom). The 
straight color differences are well within ±1 magnitude (dotted lines). Ta- 
ble[2|fits the residual dependence on IRAS [12]-[25] color. The MSX color 
does not depend on total IRAS 12 /im flux density (neither is corrected for 
zero-magnitude offset). Results are similar for IRAS 25 and 60 /im flux den- 
sities (not shown). For very red sources, IRAS [12]-[25] > 1, only the MSX 
[C]-[E] color is closely related to the IRAS [12]-[25] color. 



(Fig. |2]). Color region VIII contains a complex mixture of 
sources. For the remainder of this paper we will concentrate 
on the oxygen-rich CSE-sequence I-V, completely ignore the 
88 sources in color region VIII, and only occasionally com- 
ment on the carbon-rich sources in color regions Via, VIb and 
VII. Later we will also introduce MSX color regions with non- 
capitalized Roman numerals (/, ii, iiia, iiib, iv and v) (to keep 
them recognizable as distinct MSX, rather than IRAS color 
regions), but where the same type of sources can be found in 
the different 2CDs (Sect.@). 

In our analysis we have chosen not to apply any zero- 
magnitude or color correc tions. One of the streng t hs of the 
IRAS 2CD as defined by Ivan der Veen & Habingl (119881) is 
that it also did not have any of these correc tions applied. The 
analysis of Ivan der Veen & Habing dl988l) was therefore in- 
dependent of different choices of extinction curves and spec- 
tral ener gy d istributions (SEDs) of individual objects. Only 
in Sect. 13.21 we will temporarily divert and apply a zero- 
magnitude correction - a simple axis shift - to compare our 
data with other data sets. The zero-magnitude correction is 
applied by, e.g., adding 2.5 x log 1() (y§4§) = 1.2622 to the un- 



corrected [A]-[D] color and 2.5 x log 10 (^) = 0.7943 to the 
uncorrected [D]-[E] color, where 58.49^ 26.51, 18.29 and 
8.80 Jy are the zero-magnitude flux densities for the MSX 
bands A, C, D and E, respectively. 

Finally, although we made an effort to exclude data that is 
"confused" (i.e., where the IRAS beam would include emis- 
sion from several objects as compared to emission from a sin- 
gle object in the smaller MSX beam), it is quite possible that 
a few sources are confused or mis-identified. However, from 
the tail of Fig.Q]we estimate the number of mis-identifications 
or confused sources is less than 10%. Since we are making a 
statistical comparison of the colors, we expect no major dif- 
ference in our conclusions when a few misidentified objects 
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TABLE 1 

Least square fits" of IRAS [12]-[25] versus MSX colors 



MSX 


slope 


offset 


RMS 


MSX 


slope 


offset 


RMS 


[A]-[C] 
[A]-[D] 
[A]-[E] 


0.468 
0.852 
1.227 


0.230 
0.133 
0.140 


0.224 
0.283 
0.296 


[CHP] 

[C] -[E] 

[D] -[E] 


0.375 
0.765 
0.384 


-0.099 
-0.100 
0.003 


0.154 
0.177 
0.205 



"; MSX color = slope ■ (IRAS [12]-[25]) + offset 



remain in the sample. 

3. FROM IRAS COLORS TOWARD MSX COLORS 
3.1. The IRAS CSE-sequence in MSX Colors 

The CSE-sequence of oxygen-rich objects in the IRAS 2CD 
(see Fig. [2]) is an exponential function of the IRAS [12]-[25] 
color resulting in an IRAS [25]-[60] color. With respect to 
the [12]-[25] color region domain the CSE-sequence over- 
laps in regions Illb and IV (between 0.4 < [12]- [25] < 1.0) 
and in regions IV and V (between 1.0 < [12] -[25] < 1.4). 
However, this overlap is much less than in the [25]-[60] 
color region domain, where regions I-IIIb overlap between 
—2.4 < [25] - [60] < -1 .6. It is therefore convenient to regard 
the CSE-sequence in terms of IRAS [12]-[25] color. Using 
the [12]-[25] colors at this stage does not distinguish between 
carbon-rich objects in regions Via, VIb and VII and oxygen- 
rich objects along the CSE-seq uenc e for [12]-[25] < 0.4; we 
will discuss this further in Sect. 13.21 

Strictly speaking, in terms of the center wavelength of 
the wavebands, the IRAS [12]-[25] color would correspond 
roughly to the MSX [C]-[E] color, with no corresponding 
color in the Spitzer (GLIMPSE+MIPS) n or in the ISOGAL 
surveys. However, as readily explained in lOrtiz et al.l d2005). 
the MSX A filter covers the 9.7 /im silicate feature, whereas 
the 11.3 /im SiC feature would appear in the MSX C filter; the 
IRAS 12 /im band filter includes both features 4 . As the sili- 
cate feature is typically the most prominent in the IRAS LRS, 
the IRAS [12]-[25] color may also, or better, correspond to 
the MSX [A]-[E] color. Figures & ED and Table Q] com- 
pare the IRAS colors to the MSX colors that could be derived 
from the cross-identifications. This comparison shows that 
the IRAS [12]-[25] color can generally be converted to MSX 
colors, most notably in the MSX [A]-[E] and [C]-[E] colors, 
but also with a reasonable color range in the MSX [A]-[D] 
color (Table |2). We note that there is no dependence on total 
flux density at any waveband (Fig. [4]i. Except for the [C]-[E] 
color, the dispersion and deviation from the trend is in general 
larger for very red sources, where [12]-[25] > 1, but overall 
the agreement is well within a third of a magnitude. We note 
that most of these sources are intrinsically variable evolved 
stars for which a difference in color due to non-simultaneous 
observations within a mission, or observing different stellar 
phases between the IRAS and MSX missions is reflected in 
this dispersion dRobitaille et alJl2007l) . No conversion can be 
made for the IRAS [25]-[60] color. This is not surprising as 
all MSX colors are determined for the SED shortward of 25 
/im, and are largely insensitive to the SED and IRAS color 
longward of 25 /im wavelength. 

For the oxygen-rich sources along the CSE-sequence, this 
close correspondence to the IRAS [12]-[25] color to three 
MSX colors is fortunate, as the CSE-sequence thus can be 

4 The Spitzer InfraRed Array Camera (IRAC) 8.0 fim filter excludes both, 
which could complicate recognizing the CSE-sequence using IRAC data. 



TABLE 2 

Least square fits of IRAS/MSX color differences 



IRAS - MSX colors slope offset RMS 



(IRAS [12]- 


-[25]) - 


(MSX [A]- 


-[D]) 


0.148 


-0.133 


0.283 


(IRAS [12]- 


-[25]) - 


(MSX [A]- 


-[E]) 


-0.227 


-0.140 


0.296 


(IRAS [12]- 


-[25]) - 


(MSX [C]- 


-IE]) 


0.235 


0.100 


0.177 



represented by the IRAS [12]-[25] color, and also by the 
MSX [A]-[E], [C]-[E] and [A]-[D] colors. We will thus as- 
sume that the small inaccuracies of converting IRAS colors 
using the linear relations in Table Q] remain within the noise 
of the whole sample of derived MSX colors. Then the IRAS 
CSE-sequence can be described as MSX CSE-sequences by 
substituting the I RAS [12]-[25] color with two d ifferent MSX 
colors in Eq. 2 of lvan der Veen & Habind (fl988) for the range 
of valid IRAS [12]-[25] colors (-1.5 < [12]-[25] < +1.3). It 
follows that the IRAS [25]-[60] color can be eliminated by 
equating the right hand sides of the two equations obtained 
above. Rearrangement of the resulting eq uatio n yields the 
MSX CSE-sequences given in Tablel3l(Sect.l4~TTi. 

3.2. Carbon-rich or oxygen-rich? The MSX "OC-index" 

In a statistical sense IRAS colors can be used to sepa- 
rate the carbon-rich sources from the oxygen-rich sources, as 
the carbon-rich CSEs typically reside in IRAS regions Via, 
VIb and VII. However, when only MSX colors are avail- 
able the separation of carbon-rich from oxygen-rich envelopes 
becomes more of a problem. Although other information 
is sometimes available upon which a clearer distinction be- 
tween carbon-rich and oxygen-rich CSEs can be made, (e.g. 
J K o r L band 2MAS S or D ENIS flux d ensities, iSuh etalJ 
d200lb ; iLumsden et ail d2002l) ; ISuhl d2007h ). the challenge is 
to make such a distinct ion using only the data available within 
the MSX c atalog (e.g.. ILumsden et al.ll2002b ISevensterJl2002l; 
lOrtiz e t al. 2005|). In our view the most promising case is 
the data in Fig. 3 of lOrtiz et al.l ( 120051) wh ich implicit l y sug- 
gest that a first separation can be made. Ortiz etail (2005) 
show that discriminating between C and O rich stars is bet- 
ter done using K- [A] versus [A] - [D] colors (their Fig. 2), 
but here we only use MSX colors and turn to their analysis 
of the zero-magnitude corrected MSX [D] - \E] versu s MSX 
[A] - [D] color (their Fig 3). The s amples in Ortiz eTall (120051) 
are well defined as the data from lLoup et al.l (119931) that they 
use are identified as carbon-rich stars and indeed typically, 
but not uniquely, pop ulate IRAS color regions Via, VIb and 
VII dLoup et al.lll993l their Fig. 3). Similarly, but vice versa 
for oxygen-rich stars in regions I-V, this reminds one that the 
IRAS boundary between carbon-rich and oxygen-rich color 
regions is not strict. 

In Fig. 3 of lOrtiz et all (120051) the carbon-rich stars (filled 
squares) are mostly separated from the oxygen-rich OH/IR 
objects by the line [A] - [D] = -2 ■ ([D] - [£]) + 2.5 using zero- 
magnitude corrected colors. This separator is the dotted line 
in our Fig. [5] and labeled "C vs O" 5 . Is it possible to use this 
empirical line to generally make the division between carbon- 
rich and oxygen-rich objects when the composition of the en- 
velope is not known a-priori, e.g., such as in samples used in 

5 As an interesting side note, thes e are colors involving the MSX D filter, 
whereas in the previous section 13. It we mentioned that the "carbon feature", 
the SiC line, appears in the MSX C filter. Plausible components with carbon- 
bearing molecules which could affect the fl ux density in the D filter are those 
of HCN and C 2 H 2 at ~ 14 fim (see e.g.,[Blommaert et aU2005D . 
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Fig. 5. — Zero -magnitude correct ed M SX \D] - [£J vers us [A] - [D] 
2CD following loi-tizetal J d2005h (and lOjha etal]|2007l) for the fi- 
nal selection of 2455 sources after removing those in IRAS color 
region VIII (88 sources). For 544 sources data could not be plot- 
ted in this diagram (see text). Filled squares are sources from 
IRAS color regions Via, VIb and VII, the supposedly carbon-rich 
sources; open circles are supposedly oxygen-rich sources from IRAS 
color regions I-V. The dotted line labeled "C vs O" ([A] - [D] = 
—2 ■ ([D] - [ EJ) + 2.5) is an empirical line derived from Fig. 3 in 
I Ortiz et aT] (2005). It suggests a division between carbon-rich and 
oxygen-rich sources in their data (see text a nd Fig. H}. Also drawn 
(dashed) are the lines from lOrtiz et alj |2005) and Oiha et al.l d2007t) 
that depict a gap between stars and planetary nebulae or transition 
objects in their data; the gap is also apparent here. The (solid) black 
body temperature line ranges from 10,000 K at the origin to about 
14 K at [D] — [E] = 2 (note that the black body temperature scale 
in lOrtiz et alj (2005) is in error [R. Ortiz, priv. comm.]), crossing 
the "C vs O", "Ojha07"and "Ortiz05" line at T ~ 400, 255 and 230 
K respectively. These lines are roughly perpendicular to the black 
body temperature line and thus may be interpreted as temperature 
separators. 

lOjha et alJ (120071) and here? 

In Fig. [5] the sources of our sample in IRAS color regions 
Via, VIb and VII (filled squares) and the other regions (open 
circles) excluding region VIII (88 sources) seem to overlap 
in the densely populated part of the figure around the "C vs 
O" line. Note that only 1911 of the original 2455 sources 
could be plotted in this figure, as 544 sources are missing at 
least one good quality A, D, or E flux measurement. Fig- 
ure|6]separates the two samples to show that supposed carbon- 
rich sources from IRAS region Via, VIb and VII do appear in 
the oxygen-rich (top right) side of the "C vs O" dotted line. 
Supposed oxygen-rich sources also appear on the carbon-rich 
(bottom left) side. The IRAS color region s VIb and VII as 
described by Ivan der Veen & Hab ing ( 1988) may still contain 
a large fraction of oxygen-rich st ars based on the p ercentage 
of LRS classes 21-29 and 14-16 dOlnon et al j 19861) reported 
in the samples ( Ivan der Veen & Hab ing 1988). There is also 
some fraction of supposed oxygen-rich stars in regions II- 
IV which have LRS classes con sistent with carbo n-rich stars. 
This may also be seen, e.g., in ISuh et all ( 12001b where evo- 
lutionary models of carbon-rich stars do traverse oxygen-rich 
IRAS regions Ilia and Mb, and models for oxygen-rich stars 
border region VII. 
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FIG. 6. — Zero-magnitude corrected MSX [D]-[E] versus [A]-[D] 2CD as 
Fig.|5] but now separated for IRAS color regions Via, VIb and VII ("Carbon- 
rich"; 585 sources for which an "OC-index" could be calculated out of 966 
sources in the carbon-rich IRAS regions Via, VIb and VII) versus IRAS color 
regions I-V ("Oxygen-rich"; 1326 of 1489 sources). The (dotted) line is the 
empirical "C vs O" line from Fig. \5\ Selecting sources to the lower left of 
the "C vs O" line, i.e., selec ting sources for which the "OC-index" (Eq.[T} is 
less than zero (see text; Sect. 13.2k helps to increase the fraction of carbon-rich 
sources taken from the sample. However, selecting sources to the upper right, 
with a positive "OC-index", very effectively deselects oxygen-rich sources 
with T > 400 K (see FigTTJ. 




IRAS [121 f25l 



FIG. 7. — IRAS 2CDs of the 1911 sources separated according to the OC- 
index ("C vs O" line) for potentially carbon-rich objects (left; 886 sources) 
and oxygen-rich objects (right; 1025 sources). In general the division using 
the OC-index puts the sources in the expected IRAS regions, but the OC- 
index is not a good divider for sources hotter than 400 K (where the "C vs 
O" line crosses the black body temperature line in Fig. [3} as no oxygen-rich 
objects end up in IRAS regions I and II. 

Not surprisingly this means that the IRAS region division 
between oxygen-rich and carbon-rich stars is not perfect. We 
therefore investigated whether a better separation between 
carbon-rich and oxygen-rich objects can be made using the 
MSX [D]-[E] versus [A]-[D] colors dividing at the "C vs O" 
line, with the understanding that the MSX separation may also 
not be perfect. We define the MSX "OC-index" for sources 
with good quality MSX A, D and E flux densities: 

OC-index = ([A]-[D]) + 2 -([D] -[£]) + C, (D 

with £ = -2.5 for zero-magnitude corrected MSX colors and 
C = +0.35 for uncorrected MSX colors. The OC-index is thus 
negative if the source is on the supposedly carbon-rich side of 
the "C vs O" line and positive on the supposedly oxygen-rich 
side. 

Figure [7] shows the IRAS 2CD for the sample, divided ac- 
cording to the OC-index separately. Although the OC-index 
seems to generally divide the carbon-rich from the oxygen- 
rich objects, it mistakenly identifies potential oxygen-rich 
sources in IRAS regions I and II as carbon rich. It does se- 
lect the cooler and thicker CSEs in IRAS color regions Illa-V, 
but since we are interested in the full CSE-sequence, (includ- 
ing IRAS color regions I and II), using this division does not 
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suit our needs in this context. We conclude that regardless 
of the potential to increase the fraction of carbon-rich sources 
in samples of hot objects (T > 400 K, left of the border be- 
tween region VII and Ilia), the OC-index is more a measure 
of temperature than CSE chemistry (Pig. [7j. 

The main interest in the present study is in the oxygen-rich 
sources along the CSE-sequence. As none of the MSX colors 
relates closely to the IRAS [25]-[60] color, it is expected that 
the carbon-rich and oxygen-rich sources overlap in MSX col- 
ors where the uncorrected IRAS [12]-[25] color is less than 
~ +0.4, the farthest extent of IRAS color region VIb. To elim- 
inate the pollution of an oxygen-rich sample along the CSE- 
sequence with likely carbon-rich sources in the MSX color 
samples, we therefore removed from our sample all sources 
from the carbon-rich IRAS color regions Via, VIb and VII 
(and VIII). However, we did not remove the 388 sources with 
OC-index < from the sample as this would eliminate all hot 
(T > 400 K) oxygen-rich objects as well. Thus a reduced sam- 
ple of 1489 sources with good MSX colors in the oxygen-rich 
IRAS regions was used in the remainder of this paper. Note 
that we chose to eliminate the likely carbon-rich sources. An 
option to introduce different symbol colors (see Fig. [8) would 
not add to clarity in the figures later on, whereas opting for 
symbol colors according to their IRAS [12]-[25] color, and 
thus mislabeling them as oxygen-rich sources, would just plot 
them on top of the distribution of oxygen-rich sources with 
the same symbol color. That is, leaving them out from the 
comparison of infrared colors of oxygen-rich sources seems 
proper and does not alter our conclusions. 

4. MSX TWO-COLOR DIAGRAMS 

Using the sample of 1489 generally oxygen-rich sources, 
MSX 2CDs were constructed. In an attempt to describe MSX 
regions that would correspond to IRAS regions, it is impor- 
tant to determine MSX 2CDs in which the different IRAS re- 
gions clearly separate one from another. As with the IRAS 
regions, some overlap of sources is expected and unavoid- 
able with these statistical approaches. All possible combi- 
nations of MSX 2CDs were investigated. However , as can 
be deduced from the IRAS-MSX relations in Sect. 13.11 and 
Fig.lU the MSX colors [A]-[D], [A]-[E], and [C]-[E], which 
have a slope larger than 0.5, are most useful. Furthermore, a 
color using the E-filter is the most distin ctive as it gives the 
largest color range (but note Sect. 14.1. 11 1. Also, because all 
MSX colors represent the IRAS [12]-[25] color, a more or 
less linear relation between any of these two MSX colors may 
be expected albeit with different absolute color ranges. In 
Fig.[8]these MSX 2CDs are shown with different source col- 
ors, depending in which IRAS color region the MSX source 
resides (see top right panel). In general the trend is along the 
CSE-sequence which is no surprise as the CSE-sequence, the 
black body temperature curve and the IRAS [12]-[25] color 
all align with the sequence of IRAS color regions I-V 

4. 1 . MSX color regions and the CSE-sequence 

Figure |8] shows that sources from a particular IRAS region 
(a specific color) generally populate a confined region in the 
MSX 2CD diagrams. We therefore determined MSX 2CD 
color regions for sources with IRAS colors of color region 
I, II, etc. We will name the MSX region with sources from 
IRAS color region I, MSX color region i and similarly, e.g., 
sources from IRAS color region Mb, MSX color region iiib, 
etc. Using our limited sample of 1489 sources total, we empir- 



ically determined approximate MSX color region boundaries 
as follows (Fig. [8): 
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+0.4 < [A]-[D] < +5.0 and 

+0.1 < [C]-[E] < +0.55 or 

MSX Uib . < 1.4-([A]-[E])-([A]-[D])<+0.8 and 

+0.4 < [A]-[E] < +1.1 or 

+0.1 < [C]-[E] < +0.55 

+0.8 < 1.4-([A]-[E])-([A]-[D])<+1.3 and 
+0.4 < [A]-[E] < +7.0 or 



MSX IV 



MSX V 



+0.55 < [C]-[E] < +1.0 

1.4-([A]-[E])-([A]-[D])>+1.3 and 
+1.1 < [A]-[E] < +7.0 or 

+1.0 < [C]-[E] < +5.0 



Two regions in the [A]-[E] versus [A]-[D] diagram (up- 
per left panel in Fig. [8) are not accounted for in the above 
determination due to the lack of sources. The first region is 
for 1.1 < [A]-[E] < 7.0 and with 1.4-([A]-[E]) - ([A]-[D]) 
< +0.8, where sources from either IRAS region Mb and IV 
can be expected. Sources from either IRAS region IV and 
V can be expected between +0.4 < [A]-[E] < 1.1 and with 
1.4-([A]-[E])-([A]-[D])>+1.3. 

We also included the beginning, zero and end points of the 
CSE-sequence in the MSX 2CDs. These are the locations of 
the crosses in the MSX 2CDs (Fig.© where IRAS [12]-[25] 
equals to -1.5 (labeled C), zero (labeled S) and +1.3 (labeled 
E). The full CSE-sequence is a straight line between these 
points and can be obtained by elimin ating the [251— [60"! color 
from t he CSE-sequence equation in Ivan der Veen & H abing 
(1988, their Eq. 2), or simply by equating the expressions 
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Fig. 8. — Top right: IRAS 2CD with color coded regions along the CSE-sequence. These color codes are used for the color of each source in 
the other panels, in the MSX 2CDs. For example, sources with IRAS colors consistent with being located in IRAS color region Illb (yellow) 
appear in yellow in the MSX 2CDs. The black body temperature curve temperature labels are omitted at the ticks for clarity; they are the same 
as in Fig.|2] i.e., 10,000 K, 2000, 1000, 500, 300, 250, 200 and 160 K from left to right respectively. Other panels: MSX 2CDs of the 1489 
sources in IRAS regions I-V. The color of the source depicts in which IRAS region it resides (see top right panel). Small cyan (lighter blue) 
points are sources with good MSX colors from the entire MSX catalog. The IRAS matches discussed here are a subset of mainly brighter 
sources which are in the Galactic plane. The CSE-sequence is the line connecting the crosses labeled C, S and E (the line itself is not drawn, 
for clarity), plotted at [12]— [25] = -1.5, and +1.3 respectively. The solid line is the black body temperature curve starting at 10,000 K in the 
lower left, interrupted between 230 K and 200 K for reference, and ending at 130 K. The fat-dot ticks are at the same temperatures as in the 
IRAS 2CD. The dashed lines in the top left panel depict the lines separating sources from IRAS region Illb from those from IRAS region IV 
and sources from region IV from those from IRAS region V - see text. 



for [12]-[25] from Table Q] for different MSX colors over 
the [12]-[25] range in MSX colors where the CSE-sequence 
is defined. The expressions and color ranges for the CSE- 
sequence in MSX colors adapted here are given in Table [3] 
It can be seen that indeed the CSE-sequence starts in MSX 
color region i, has its zero point close to the border between 
regions iiia and iiib 6 , and ends in region iv. It does not en- 

6 In Fig. |s] the yellow symbols of region iiib overlap and unfortunately 
overplot the blue symbols of region iiia, making it unclear that the blue sym- 
bols continue under the yellow symbols. However, note that the yellow sym- 



ter the region (v) with the reddest sources, which are mostly 
(proto-)planetary nebulae. 

To obtain a subsample of sources that is characterized by 
colors of a certain IRAS color region, the selection criteria 
for the same (lower case) MSX color region can be directly 
applied if the parent sample contains uncorrected MSX col- 
ors. For example, a subsample of sources with CSEs re- 
sembling the IRAS Ilia region colors (blue in Fig. [8]l can 
be drawn directly from the MSX catalog by applying one 

bols do not extend much past the negative (lower left) side of [A]-[E] = 0. 
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FIG. 9. — Resulting IRAS 2CD when the selection for MSX color regions 
is applied. IRAS color regions I, II and Ilia can readily be selected using a 
single MSX color, [A]-[D], [A]-[E] or [C]-[E]. Because the MSX [C]-[E] 
color corresponds best to the IRAS [12]— [25] color, IRAS color regions Illb, 
IV and V can also be selected using the [C]-[E] color. If the [C]— [E] color is 
unavailable, due to a missing or unreliable C filter flux density, these regions 
can only be selected using a combination of the [A]— [D] and [A]— [E] colors. 

or more of the three equations for MSX region iiia. On the 
other hand, to obtain such a subsample from a parent sam- 
ple that contains zero-magnitude corrected MSX colors, one 
woul d hav e to subtract the zero-flux magnitude fractions from 
Sect. 12.21 from the MSX source colors first 7 . As a reminder, 
we emphasize that carbon-rich sources in IRAS regions Via, 
VIb and VII will overlap with MSX color regions i, ii and 
iiia, with a small tail in iiib-v. Using only MSX colors it 
is impossible to make an accurate estimate of the fraction of 
contamination with carbon-rich sources in a sample. How- 
ever, to get an indication, in Table [4] we show the number of 
sources from the cross-identified IRAS/MSX data base (2543 
sources; Sect. |2]i in each IRAS color region, when sources 
are selected using any of the MSX selection rules above. It 
is apparent from this table that selecting a sample using the 
MSX color relations in general selects the sources with simi- 
lar IRAS colors (e.g., MSX iiia and IRAS Ilia), but that also 
sources from nearby IRAS regions may fulfill the MSX se- 
lection rule. This, in a statistical approach, is not worrisome; 
e.g., using the MSX [A]-[D] color to select MSX iiia sources 
will yield 82% (816/995) of IRAS Ilia sources, provided that 
the sample is known to be mostly oxygen-rich along the CSE- 
sequence. Because of the degeneracy in the vertical axis of 
the IRAS 2CD when only MSX colors (or IRAS [12]-[25] 
colors) are available, the carbon-rich objects in IRAS regions 
Via, VIb and VII also fulfill the MSX selection rules. The 
approximate number of sources in these IRAS regions (ac- 
cording to increasing [12]-[25] color) is also given in Ta- 
ble [4] per MSX selection rule. If the degeneracy between 
oxygen-rich and carbon-rich cannot be resolved, e.g. with 
the IRAS [25]-[60] color, then the pure fraction of IRAS re- 

7 Because the MSX A and D waveband filters are similar to the ISOGAL 
7 and 15 (im waveband filters, selection on the zero-magnitude corrected 
ISOGAL [7]— [15] color ranges may appear to be si milar to selection on th e 
zero-magnitude color corrected MSX [A]— [D] color (Messineo et al. 2004). 



TABLE 3 

CSE-SEQUENCE IN MSX TWO-COLOR DIAGRAMS (FlG.[8]l 

MSX CSE-sequence relations Valid MSX CSE color ranges 

[A]-[D] = 0.694 ■ ([A]-[E]) + 0.036 -1.145 < [A]-[D] < +1.241 

[C]-[E] = 0.623 ■ ([A]-[E]) - 0.187 -1.701 < [A]-[E] < +1.735 

[C]-[E] = 0.898 ■ ([A]-[D]) - 0.219 -1.248 < [C]-[E] < +0.895 



gion Ilia sources of the MSX iiia sources using the [A]-[D] 
color selection degrades to 50% (816/1624). Nevertheless, 
selecting on e.g. MSX iiia will thus indeed mostly yield 
IRAS Ilia sources, partly because some of the IRAS Via, 
VIb and VII sources may be actually oxygen-rich since the 
IRAS boundaries also do not exactly separate the two popu- 
lations. Because the density of carbon-rich sources is highest 
for [12]-[25] < 0, the highest contamination (cq. overlap) is 
in MSX color regions i, ii and iiia. For these regions other se- 
lection crit eria to increase the o xygen-rich fraction are advis- 
able (e.g. Messi neo et alJ (120021) . Sect.|5]l. Again, the oxygen- 
rich and carbon-rich objects overlap around the dividing re- 
gions in the IRAS 2CD, in particular near IRAS region VII, 
so the numbers in Table [4] are indicative values only for any 
sample other than this IRAS/MSX cross correlated data set. 
For the particular example above, selecting MSX iiia sources 
using the [A]-[D] color, the fraction of oxygen-rich objects 
is 61% (995/1624). That is, for such a sample the remainder, 
39%, may be carbon-rich and thus will not harbor SiO masers 
to begin with (Sect. [5]). 

4.1.1. Possible error in MSX PSC 2.3 ? 

In our investig ation we noticed possible error in the MSX 
PSC version 2.3 dEgan et al.ll2003l) . In the MSX [A]-[E] ver- 
sus [A]-[D] 2CD (the top left panel in Fig. [HJ there seems 
to be a different distribution of sources which is shifted to 
the right and away from the CSE-sequence with respect to the 
main-populated area when all MSX sources with good quality 
photometry are used 8 . We could not account for this different 
distribution except by positing a somewhat different, redder 
population. In order to help confirm this possibility, we made 
the same 2 CD diagram with data retrieved from MSX PSC 
version 1.2 (Eg an et al.lll999t) . When compared to the PSC 
version 2.3 2CD there seems to be a more continuous distri- 
bution of redder (and weaker) objects in the diagram made 
from version 1.2 data (not shown). 

From both versions of the PSC we selected objects with 
quality 3 or 4 in A, D, and E and constructed each color from 
the catalog data. Typically the sources weaker in flux density 
in the MSX A filter are affected, with a shift in [A]-[E] of 
~ 0.5-1 magnitude redward. We understand that the differ- 
ence between PSC versions 1.2 and 2.3 is that the photometry 
has been improved, and weaker sources have been included. 
However, we would not have expected the "gap" in the distri- 
bution of colors in the 2.3 data (top left panel in Fig. [8J that 
does not appear in the 1 .2 data. The reason for this is currently 
unknown, but seems to originate from MSX E filter (21.3 fj,m 
band) fluxes that are a factor of about two too high for most 
of the weaker A filter (8.28 /im band) sources 9 . 

8 This may also be seen in the [A]-[D] versus [C]-[E] 2CD in the bottom 
right panel, where many sources show up above the main-populated area. 

9 A referee has explained the probable cause to be a combination of un- 
compensated variable response during the mission and the different tech- 
niques used for the extraction of the photometry in the different versions. 
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TABLE 4 

Correlatio n ma trix for sources for which the MSX color could be calculated and fall within the ranges of the equations 

GIVEN IN SECT. |4.1| ("T0TAL"). THE NUMBERS ARE THE SOURCE COUNTS IN EACH IRAS COLOR REGION FOR EACH OF THE MSX COLOR SELECTION 
RULES, AND INDICATE THE POTENTIAL CONTAMINATION LEVEL OF OTHER SOURCES THAN TARGETED USING THE MSX RULES. 
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We have investigated whether our analysis also suffers from 
possible inaccuracies or effects in the weaker 8.28 /im ob- 
jects. Fortunately the cross-identified objects, that are also 
detected by the less sensitive IRAS 12 [im detector, belong 
to the brighter MSX sample for which almost no sources are 
expected to be affected. This can also been seen from the lack 
of color-coded source symbols shifted to the right of the main 
distribution in the [A]-[E] versus [A]-[D] diagram. Note that 
the black dots of MSX color region v (or IRAS color region 
V, i.e., mostly planetary nebula) are shifted redward, because 
they are intrinsically redder than the CSE-sequence sources. 

5. 86 GHZ SIO MASERS AND MSX COLOR 

Previously the IRAS 2CD has been very successful in se- 
lecting objects to target: for example, for 1612 MHz maser 
emission (e.g. Jte Lintel Hekkert et al.ll991l) . Similarly, IRAS 
colors have been used to select targets for 86 GHz (J = 2 — > 
1 , v = 1) and 43 GHz (J = 1 -> 0, v = 1 and v = 2) SiO maser s 
(e.g..lHaikalal[T990t lHaikala et al.l[l994l iNakadaet al.l[l993l) . 
Jiang et al. ( 1995J)have made a careful analysis of the 43 GHz 
SiO maser detection rate as function of IRA S color and IRAS 
color r egion from the data collected by llzumiura et al.1 (119941 
Il995allbh in the Galactic bulge and disk. The Izumiura data 
were selected based on color and variability index with likely 
carbon-rich stars from IRAS regions Via, VIb and VII re- 
moved. 

In addition to an overall constant detection rate of SiO 
masers of about 61±3% as function of [12]-[25] color, 
I Jiang et all ( 1 19951) suggest that the detection rate in the Galac- 
tic bulge is higher than in the disk. Since in both the Galactic 
bulge and disk the IRAS survey suffered from confusion, the 
Jiang et al.l d 19951) study may have suffered from selection ef- 



fects. It will be interesting to repeat this analysis based on 
MSX colors alone. 

Another approach was used by Messineo et al. (2002)) who 
bypassed the IRAS color selection and used the combined 
ISOGAL-DENIS data to empirically determine a compli- 
cated, extinction-corrected, color-magnitude-based selection 
from MSX [A]-[D] color and also on the MSX variability in- 
dex in the Gal a ctic b ulge region. A cross-identification of 
iMessineo et all (120021) ' s target sources with the IRAS PSC 
(10" radius as in Sect. [3J results in only three sources with 
good photometry in all three IRAS bands, and only 8 1 with 
good photometry in [12]-[25] color. A selection based on 
IRAS color in this case would not have yielded an interesting 
sample, whereas the simple MSX tool presented here would 
have been extremely helpful in ease of use compared to their 
selection. 

To assess their selecti o n in t erms of the [12]-[25] color 
analysis by iJiang et al.l d 19951) . we plot the data from 
IMessineo et all (120021) in the MSX 2CDs in Fig.QHand (non)- 
detection histograms in Fig. [TT] Most sources selected by 
IMessineo et al.l (120021) are located in MSX region iiia which 
was the IRAS region they were targeting. Their sample 
likely includes m ore carbon-rich sources t han the sample of 
Jiang et al.l i 19951) as Messine o et al.l (120021) could not remove 
those s ources using the IRAS [25]-[60] color. 

The Messine o et al.l (120021) 86 GHz SiO maser detection 
rate is consistent with an overall constant 65-66% as a func- 
tion of either [A]-[D], [A]-[E] or [C]-[E] color, i.e., with 
(unavailable) f 121 — [251 color (Fig. [TTl i. Remarkably, this is 
the same result as I Jiang et all ( 1 19951) for the 43 GHz masers. 
For M SX region iiia, where the majority of the lMessineo et al.l 
d2002T) sample lies, the IRAS/MSX cross-identified sample 
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FIG 10.— MSX 2CDs for 86 GHz SiO maser targets from Messi neo et all 
(2002). Crosses are detections, open circles are non-detections. Most sources 
selected by [Messineo et al. ( 20021 are assigned to MSX color region iiia, 
comparable to the IRAS c olor region (Ilia) w here for the Galactic bulge many 
SiO masers are expected I Jiane et al. 1995). Dotted ho rizon tal and vertical 
lines show the MSX color region boundaries from Sect. 14.11 The other lines 
are as in Fig. [8] 

(Sect. [2]i indicates that the carbon-star contamination is about 
23-39%. Therefore, when selecting MSX iiia sources, us- 
ing [A]-[D] colors (similar to the Messineo sample), we find 
6 1 % to be in the oxygen-rich region of the IRAS 2CD (see Ta- 
ble Assuming that carbon-rich sources have all their oxy- 
gen locked up in carbon-monoxide molecules and thus never 
harbor SiO masers, this suggests that 1) almost all oxygen- 
rich sources selected using iiia would show detectable 86 
GHz SiO maser e mission, and 2) the com plicated selection on 
[A]-[D] color by Messin eo et alj (|2002), which includes flux 
and variability information, was likely not much better than 
simply selecting on MSX color as they both yield very sim- 
ilar detection rates. As individual histogram bins at the blue 
(or red) end of the histograms in Fig. QT]typically indicate an 
80% (or 50%) detection rate, there may be a small trend of 
decreasing SiO maser detection rate for redder sources. How- 
ever, the range in MSX color is too small, and the uncertain- 
ties too large, to make a definite statement ab out this. 

Whe n compared to the MSX PSC 2.3, iMessineo et all 
(2002)'s list of targets resulted in 204 cross-identifications 
among the 271 SiO maser detections and 109 cross- 
identifications among the 173 targets without detected SiO 
maser emission 10 . The targets without cross-identifications 
must be targets that could only be selected using the more 
sensitive ISOGAL data - thus a sample of generally weaker 
infrared sources. As relatively more SiO maser detections are 
found in the brighter infrared MSX selection this probably 
confirms that brighter infrared sourc es are more likely to sh ow 
detectable SiO maser emission (e.g., Messi neo et al.l l2002). 

A more careful study can be done with a larger sample, us- 
ing a larger color range. Here we have shown how such a 
sample selected on MSX colors (which are directly related to 
IRAS colors in the Ivan der Veen & Hab ing (1988) 2CD) can 
be obtained quite easily compared to previous selections in 




-1 o 1 

MSX color 



FIG. 1 1 . — Detection r ates as a function of MSX color for the SiO masers 
in Messineo et al. (2002). Shown are the histogram of targeted sources and 
the shaded histogram of the detected sources. The detection rate per 0.2 mag- 
nitude bin is depicted by dots, which are connected and have error bars only 
if the number of targets in the bin is ten or more. We note that the MSX 
photometry has changed slightly between PSC (preliminary) version 1 .2 and 
version 2.3, so Messineo et al. 1 2002)'s strict [A]— [D] color selection, and 
the inclusion of ISOGAL selected sources, does not hold any more in our 
[A]-[D] color diagrams for some individual sources. 

areas where IRAS was confused. 

6. POSSIBLE APPLICATIONS AND FUTURE PLANS 

It should not be too difficult to find possible applications for 
this new tool. Almost any project using IRAS color region se- 
lection toward CSEs in regions on the sky where MSX data is 
superior to IRAS data can make use of this new tool. More- 
over, particularly in the Galactic plane, new bulge and center 
projects can be attempted that were previously impossible be- 
cause of the confusion in the IRAS data. For example, one 
should be able to find a population of (proto-)planetary nebu- 
lae hidden deeply in the nucleus of the Galaxy and study how 
their CSEs are affected by the more hostile local environment. 
Another example would be to select thousands of targets to 
survey for SiO maser emission in the Galactic plane. Such 
a survey would e.g. yield a map of Galactic rotation based 
on a few thousand point masses; gas is much more suscepti- 
ble to non-gravitational forces. Stars and CSEs in regions too 
crowded for IRAS can be much better classified using MSX 
colors in terms of IRAS colors, etc. We leave it to the readers 
to find further applications in their own areas of expertise. 

A straightforward future plan would be to compare and de- 
rive similar tools usi ng the all-sky AKARI mission survey 
(formerly ASTRO-F; iMurakami etail 120071) with respect to 
the all-sky IRAS survey, but at much higher angular reso- 
lution. The AKARI PSC is expected to become available 
sometime after the summer of 2009. We have already com- 
mented on extending this study to the combined GLIMPSE 
and MIPSGAL data. We have begun to cross-identify the 
GLIMPSE and MSX PSCs, but the MIPSGAL (24 /xm) PSC 
needed to obtain Spitzer data longward of 8 /im has been un- 
fortunately delayed beyond the original release date. 



10 Instead of the 313 total sources here, Messineo et al. 12004) found 379 
cross-identifications total within 15" radius using MSX PSC version 1.2. 
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7. SUMMARY 

Using cross-identifications between the IRAS and MSX 
source catalogs we have derived color ranges and color re- 
gions in MSX 2CDs that correspond to color regions along 
the "evolutionary track" or the CSE-sequence in the IRAS 
2CD. Oxygen-rich objects in MSX color regions ;, ii, iiia, 
iiib, iv and v can be expected to have similar CSEs and CSE 
properties as oxygen-rich objects in IRAS color regions I, II, 
Ilia, Iiib, IV and V, respectively. This tool should be useful 
in selecting sources with different CSEs, especially in regions 
where selection using IRAS data is unsatisfactory such as the 
(IRAS) confused Galactic bulge and plane regions. 

However, it is impossible to distinguish between carbon- 
rich and oxygen-rich sources based on MSX colors alone; e.g. 
IRAS or near-infrared K and/or L band data is needed for that. 
This means that when only considering MSX colors, carbon- 
rich objects and oxygen-rich objects populate the same MSX 
regions, in particular i, ii and iiia. Selecting sources with a 
negative OC-index will increase the fraction of carbon-rich 
objects, but also deselect the low (T < 400 K) temperature 
carbon-rich objects. The reverse selection is not likely: that 
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one increases the fraction of oxygen-rich objects in a sample 
by selecting sources with positive OC-index, unless only cool 
sources are considered. 

We also seem to have revealed an error in the photometry of 
the MSX E filter band. The flux density seems a factor ~ 2 too 
high, or ~ 0.5-1 magnitude too bright in [E], for weak A filter 
band sources in the MSX point source catalog version 2.3. 
This error does not seem to be evident in MSX PSC version 
1 .2, which has a smaller number of weak sources. The reason 
for this discrepancy is probably due to the band E detector and 
a different handling of the photometry. 
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